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Current-induced torques are commonly used to manipulate non-coUinear magnetization configu- 
rations. In this article we discuss current-induced torques present in a certain class of coUinear mag- 
netic systems, relating them to current-induced changes in magnetic anisotropy energy. We present 
a quantitative estimate of their characteristics in uniform strained ferromagnetic (Ga,Mn)As. 
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I. INTRODUCTION 

The interplay between transport currents and magne- 
tization dynamics continues to be a major research topic 
in ferromagnetic metal spintronicsFi The current under- 
standing of this class of phenomena has been derived 
mainly from numerous studies of spin-transfer torques 
(STTs), which arise when spin polarized currents tra- 
verse non-collinear magnetic systems. STTs can be ex- 
ploited to achieve current-induced magnetization reversal 
and current-induced domain-wall motion, both of which 
have potentially important technological applications. 

There have been comparatively few studies of the influ- 
ence of transport currents on magnetization in uniform 
ferromagnets, presumably because spin transfer torques 
vanish in these systems. Yet, as pointed out indepen- 
dently by several researchers^i^ii^, current-induced reori- 
entation of magnetization does occur in some uniform 
ferromagnets. The first experimental fingerprint of this 
phenomenon was uncovered by Chernyshov et al^ who 
demonstrated that an electric current alters magnetiza- 
tion reversal characteristics in strained (Ga,Mn)As films 
with a single magnetic domain. 

STTs can be considered to be one member of a family 
of current-induced torque (CIT) effects by which trans- 
port currents influence magnetization in ferromagnetic or 
antiferromagnetio^ systems. The aim of this paper is to 
contribute to the theoretical analysis of current-induced 
torques in uniformly magnetized ferromagnets. 

In Sec. II we study the effect responsible for this type 
of torque, which we refer to as the ferromagnetic in- 
verse spin-galvanic effect i^ii In non-magnetic conductors 
the inverse spin-galvanic effect (ISGE) refers to current- 
induced spin density. Since a non-zero spin-density al- 
ready appears in the equilibrium state of a ferromagnet, 
the ferromagnetic inverse spin-galvanic effect has a dis- 
tinct experimental signature. Specifically, we find that 
in gyrotropic ferromagnets the magnetization direction is 
altered by a steady-state transport current. At a concep- 
tual level, we associate this reorientation with a change 
in magnetic anisotropy in the presence of a transport 
current. An important implication of this connection 
is that the magnetic anisotropy energy in the transport 
steady state of a ferromagnet which exhibits the ISGE 
is not invariant under magnetization reversal, essentially 



because the applied current breaks time reversal invari- 
ance. At a practical level, we provide a concise ana- 
lytical expression for the current-induced change in the 
magnetic anisotropy. This expression is suitable for eval- 
uation from first principles because it requires the knowl- 
edge of only the band structure of the ferromagnet and 
the lifetime of the Bloch states. At a technical level, 
our theory allows for the spatial inhomogeneities that 
inevitably occur in the magnitude of the ferromagnet 's 
exchange field at atomic lenghtscales. 



In Sec. Ill we carry out quantitative calculations for 
the ISGE of strained (Ga,Mn)As using a 4-band Kohn- 
Luttinger model. This calculation directly addresses the 
experiment by Chernyshov et aZ.— and corroborates their 
interpretation of the data. By computing the anisotropy 
field both in absence and in presence of an electric cur- 
rent, we find that in (Ga,Mn)As magnetization reversal 
may in principle be achieved solely by electric means: 
the required critical current densities are in the order of 
10^ — 10''A/cm^ and depend on the strain, Mn concentra- 
tion and hole density. Sec. IV contains a brief summary 
and presents our conclusions. 



The main conclusions of our work coincide with those 
reached by Manchon and Zhang in their independent and 
previously published work described in Ref. Q. Yet, 
our analysis highlights aspects that have not been em- 
phasized previously. First, we assert that in ferromag- 
nets with inversion symmetry, the current-induced spin- 
density vanishes to all orders in the strength of the spin- 
orbit interaction. Second, when evaluating the current- 
induced spin polarization we include a contribution from 
intcrband coherence which can become quantitatively im- 
portant in disordered ferromagnets such as (Ga,Mn)As. 
Third, we identify the current-induced transverse spin- 
density associated with the ISGE in ferromagnets as a 
consequence of a change in magnetic anisotropy in the 
presence of an electric current. We thus promote trans- 
port currents to the same status as temperature*^, gate 
voltages^iiSiiii, strai n^^i^^ and chemical processes^^, all 
of which are well-established control parameters for the 
tuning of magnetic anisotropy. 



2 



II. THEORY OF THE FERROMAGNETIC 
INVERSE SPIN-GALVANIC EFFECT 

In non-magnetic metals or semiconductors that are gy- 
rotropic, i.e. non-centrosymmetric and chiralfi^ a DC 
charge current is generically accompanied by a non-zero 
spin polarizationi^ This phenomenon is sometimes re- 
ferred to as the inverse spin galvanic effect (ISGE)ii Be- 
cause of the advent of spintronics and subsequent at- 
tempts to control spin polarization by electric means, 
even in paramagnetic materials, the ISGE has received 
widespread experimental^ and theoreticalii attention. 
The ISGE is purely a consequence of symmetry since i) 
current, which is odd under time reversal, is the dissi- 
pative response of a conductor to a DC electric-field, ii) 
spin is also odd under time reversal and therefore allowed 
as part of the dissipative response, and iii) axial vectors 
(like spin) and polar vectors (like current) are coupled 
in gyrotropic materialsJ^ The direction of the carriers' 
spin is determined by the direction of the electric field 
as well as by the axis along which inversion symmetry is 
broken 

The ISGE is sometimes viewed as a possible route to- 
ward the development of spintronics effects in paramag- 
netic materials that are as robust as effects like giant 
magnetoresistance that occur only in ferromagnetic ma- 
terials. Partly because spin-orbit interactions tend to 
be fairly weak, it appears to be difficult to make spin- 
galvanic effects in normal metals useful. In this section 
we turn the tables on this strategy by concentrating on 
the inverse spin-galvanic effect in magnetic conductors. 

In uniformly magnetized ferromagnets with inversion 
symmetry, the transport current is spin polarized because 
the conductivities of majority and minority spin chan- 
nels are different. This familiar fact is unrelated to the 
ISGE. Since spin-polarization is already present in the 
thermodynamic equilibrium state of a ferromagnet, the 
ferromagnetic ISGE is manifested not by the presence of 
a non-zero spin-density but instead by a change in mag- 
netization direction in the non-equilibrium steady-state 
which is dependent on the magnitude and direction of 
the electric field. In this paper we formulate a theory 
of the ISGE in ferromagnets by evaluating the torque 
which acts on the collective magnetization of a magnetic 
conductor due to spin-orbit interactions in the presence 
of a transport current. When the current is set to zero, 
the torque we evaluate vanishes along easy (and hard) 
magnetization directions and is normally viewed^S as a 
precessional torque due to magnetocrystalline anisotropy 
fields. These torques are in turn associated with the 
magnetization-direction dependence of the magnetocrys- 
talline anisotropy energy. At zero current, the anisotropy 
torques must change sign when the magnetization direc- 
tion is reversed because time reversal symmetry requires 
that the anisotropy energy be invariant under reversal. 
The ferromagnetic ISGE in gyrotropic crystals may be 
viewed as a change in anisotropy torque due to a trans- 
port current. Significantly, the ISGE torques are not odd 
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FIG. 1: Feynman diagram that encodes tlie transverse spin 
density induced by a current (ferromagnetic ISGE effect) 
which results in a change in the steady-state magnetization 
direction, a and h axe band labels for the quasiparticle and 
the quasihole. 



under magnetization reversal. 

The ferromagnetic ISGE is reminiscent of the magne- 
toelectric phenomena that have been extensively studied 
in multiferroic materialsj^i i.e. materials in which mag- 
netism coexists with ferroelectricity. A common char- 
acteristic of multiferroic perovskites is the presence of 
canted magnetism that stems from the Dzyaloshinskii- 
Moriya interaction. Since the direction of canting is 
determined by the symmetry of the crystal, one can 
envisage»2^ scenarios in which an electric-field-mediated 
reversal of the ferroelectric polarization causes a simul- 
taneous reversal of the canting angle or of the magne- 
tization. Another interesting property of multiferroic 
materials is the coupling between ferroelectricity and 
antiferromagnetismi^ This coupling makes it possible 
to switch the magnetization of an exchange-biased fer- 
romagnet by the application of an electric field. In spite 
of the contextual similarities, there are fundamental dif- 
ferences between the aforementioned phenomena and the 
ferromagnetic ISGE. For one thing, ferroelectricity oc- 
curs only in insulators while the ISGE occurs only in 
conductors. 

We evaluate the ferromagnetic ISGE microscopically 
within the framework of linear response theory (Fig. 

m))-- 

Ss' = xiEE\ (1) 

where 5s^ is the current-induced spin density (i S 
{x, y, z}), E is the applied electric field, and x the dissi- 
pative spin-current response function: 

Ztt ^ — ' ■ 

This linear response theory expression applies for time- 
independent uniform applied electric fields, and may be 
derived in the standard way^ by analytically continuing 
the imaginary^^ part of limi^^o Xse/'^^ where 

Xs,E = -T E Sa,6(k)<a(k)Gk,a(«t^n)Gk,b(jt^„+iw), 

(3) 

LUn = (2n -I- l)7rT is the Matsubara frequency at temper- 
ature T and lu is the frequency of the external field. In 
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Eq.([2]) sl^ij{k) and vl ^{k) are the k-dependent matrix- 
elements of the spin and velocity operators (Oa,6(k) = 
(a, k|0|6, k)) between Bloch states (|a, k)) in bands a and 
b. Note that the Bloch states are in general spinors in 
which orbital and spin degrees of freedom are entangled. 
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l/(eF — Ek.a + (— )j/2'''k,a) is the retarded (ad- 



vanced) Green's function evaluated at the Fermi energy 
ep, and Tk,a is the quasiparticle lifetime. For simplic- 
ity we have ignored disorder vertex corrections to both 
velocity and spin operators. In the numerical calcula- 
tions discussed in Section III we will in addition take the 
quasiparticle lifetime to be a phenomenological parame- 
ter which is independent of momentum and band labels. 

As we discuss below, the transverse components of 
the spin-density are directly related to the anisotropy 
field, which exerts a torque on the macrospin. On the 
same footing, the current-induced contribution to the 
transverse spin density is directly related to the current- 
induced contribution to the anisotropy field. 

For a ferromagnet with inversion symmetry xs,e = 
irrespective of spin-orbit interaction strength, for essen- 
tially the same reasons as the ISGE vanishes in normal 
conductors with inversion symmetryi^l This property can 
be verified by recognizing that in presence of inversion 
symmetry the Hamiltonian of the ferromagnet is invari- 
ant under k — k, which implies that Gk = G_kj 
Sa, h(k) = SQ^5(-k) and Wa,h(k) = -VQ,6(-k). Conse- 
quently, the right hand side of Eq. ^ vanishes after 
summing over all k. From a crystal symmetry classifica- 
tion standpoint there are 21 non-centrosymmetric crys- 
tal classes, among which three (T^, C3/1 and D^h) are 
not gyrotropic. The occurence of the ISGE is therefore 
restricted to 18 crystal classesii 

The main objective of this section is to relate the ferro- 
magnetic ISGE to a current-induced change in the mag- 
netic anisotropy field, yet before we do so it is benefi- 
cial to pave the way by reviewing the nuances of mag- 
netic anisotropy in electric equilibrium. In the absence 
of currents, magnetic anisotropy describes the depen- 
dence of the free energy of a ferromagnet on the direc- 
tion of its magnetization)^ Magnetic anisotropy origi- 
nates fromS^ magnetic dipolar interactions and spin-orbit 
interactions. The former lead to shape anisotropy in 
non-spherical samples while the latter produce magne- 
tocrystalline anisotropy by communicating the lack of ro- 
tational symmetry in the crystalline lattice to the spin de- 
grees of freedom. In practice, magnetic anisotropy reveals 
itself in dynamical processes such as ferromagnetic reso- 
nance through an anisotropy field that forces the magne- 
tization to precess unless it is along an easy or hard axis, 
i.e. along a direction in which the anisotropy energy is 
minimized or maximized. This precessional magnetiza- 
tion dynamics is properly characterized by the Landau- 
Lifshitz equation, dt^ = ^ x Hes, where Cl is the direc- 
tion of the ferromagnet 's collective dynamical variable 
(which may be chosen to be either the magnetization or 
the ferromagnetic exchange field) and H^s is an effective 
magnetic field, taken here to include reactive as well as 



riooi 




[010] 



FIG. 2: Cartoon of a magnetic thin film (shaded area). The 
exchange field A is an effective magnetic field which is parallel 
to the magnetization only when it points along easy or hard 
crystalline directions. The orientation of A can be specified 
by the polar and azimuthal angles 9 and 0. The relationship 
between the direction of A and the direction of magnetization 
is altered by an electric current in gyrotropic ferromagnets. 



dissipative processes j^^*^ The anisotropy field may then 
be defined as the contribution to the non-dissipative part 
of the effective magnetic field which survives in the ab- 
sence of true magnetic fields: 



1 dEas 
So dCl 



(4) 



where Eqs is the ground state energy of the ferromagnet 
in equilibrium (we take zero temperature throughout) 
and 5*0 is the total spin (magnetization x volume) of the 
ferromagnet. 

When we discuss (Ga,Mn)As in the following section, 
we will use spherical coordinates (Fig. [TT| in which the 
anisotropy field may be written as 



Han — Hij,(f> + HgO, 



(5) 



where and 6 are the azimuthal and the polar unit vec- 
tors, respectively. The longitudinal component of the 
anisotropy field is irrelevant because Cl y. Q, — Q. 

In order to elaborate on the microscopic theory of the 
anisotropy field in a concrete way we work within the 
spin-density-functional theory of a magnetic material, in 
which the effective Hamiltonian that describes the the- 
ory's Kohn-Sham quasiparticles can be expressed as 



H = Tiki, 



Us 



(6) 



In Eq. dB]) A = Ao(r)r2 is the exchange effective- 
magnetic- field of the ferromagnet, O is the direction of 
the exchange field, s is the quasiparticle spin operator, 
TLso captures spin-orbit interactions, and TYkin collects 
all spin-independent terms in the Kohn-Sham Hamilto- 
nian. In this work we characterize the macrostatc of a 
ferromagnet by specifying the direction of the exchange 
field. is assumed to be uniform in space but the mag- 
nitude Aq (r) of the exchange field is allowed to have spa- 
tial dependence at the atomic lengthscale.— We neglect 
dipolar interactions since they are not directly infiuenced 
by currents and can normally be cleanly separated from 
magnetocrystallinc anisotropy. 
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It follows that the zero-temperature anisotropy field is 
jiven by 



Ha 
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de 



k.a 



k,a 



(7) 



In Eq. d?]) we have used^i Eqs = I]k,a ^k,a/k,a, where 
ek,a is the energy of the Bloch state quasiparticles and 
/k,a = 9(eF — fik.a) IS the equilibrium occupation factor 
at zero temperature. Furthermore we have exploited the 
fact that 



E 



Ek.a- 



an 



(8) 



since the number of electrons in the ferromagnet is invari- 
ant under rotations of the magnetization. This implies 
a fi-dependence of the Fermi energy ,^^'^'^ which is taken 
into account in the calculations of Sec. III. 

Eq. ^ may be rewritten in a more informative man- 
ner using the Fcynman-Hcllmann theorem, which implies 
that 



Then, 



(a,k|^|a,k) 



(a,k|Ao(r)s,|a,k). (9) 



^ E^"' '^l^o(r)s|a, k)/k,a, 



So 



(10) 



where (a, k|Ao(r)s|a, k) = / drAo(r)(a, k|r)s(r|a, k). 

For the envelope-function model we use in the next sec- 
tion, the magnitude of the exchange field is a spatially 
constant Aq and the torque exerted by the anisotropy 
field is simply equal to the Aq times the transverse spin- 
density divided by the total spin of the ferromagnet. In 
ah initio calculations, the magnitude of the exchange field 
always varies substantially on an atomic scale and, as we 
have emphasized previously,— this variation must be ac- 
counted for. In this case the anisotropy field is evaluated 
by integrating the product of the exchange field magni- 
tude and transverse spin density over space. 

Eq. (jlOp may be separated into azimuthal and polar 
components: 

Jo . — 

k,a 



k,a 



If we neglect spatial variations of Ao(r), Eqs. (fTO|) and 
pT|) indicate that the torque created by the anisotropy 



field will vanish when the (spin) magnetization E(^)/ 
is parallel to the exchange field. Conversely, whenever 
the direction of magnetization is misaligned with A, the 
anisotropy field will be nonzero and will produce a torque 
on the magnetization. In transition metals spin-orbit in- 
teractions produce a misalignment between the exchange 
field and the magnetization, unless SI is pointing along 
some special crystalline direction that corresponds (by 
definition) to an easy or hard axis. A similar picture ap- 
plies to local-moment ferromagnets as well, where due to 
spin-orbit coupling the direction of the local moments is 
generally misaligned with the direction of the itinerant 
spin density. 

One of the targets of this section is to present formulae 
that are useful for researchers working on both model sys- 
tems as well as ab-initio electronic structure calculations. 
Therefore, we digress to explain that Eq. pO|) is equiva- 
lent to the alternative expressions found in ab-initio stud- 
ies. In first principles magnetic anisotropy theory^i^i 
Eq. ([9]) has been approached from a different perspec- 
tive. In such approach it is customary to choose the spin 
quantization axis along the direction of magnetization, 
so that A • s = AqSz is independent of SI. When this 
choice is made, the spin-orbit term in the Hamiltonian 
becomes explicitly f2-dependent. Consequently, 



5ek,a 



I dHso 



(12) 



The anisotropy field is then evaluated combining Eq. (fT^ 
with the force theorem^ and a full-potential electronic- 
structure calculatioui^ Of course, the final result is in- 
variant with respect to the choice of the spin quan- 
tization axis. In order to prove the equivalence of 
Eqs. (HI) and it is convenient to rewrite^ Eq. 

as de/dcf) = (90[exp(is • 2^)7^80 exp(— is • z(/))])|o and 
dejdQ ~ (9e[exp(is • 00)7Yso cxp(— is • 00)])|o. To see 
that these expressions agree with Eq. (jlip note that 
[Hso,s] = \H- Hi^in + A • s,s], that [Hkin,s] = 0, and 
that (a,k|['H,s]|a,k) (ek,a - ek,a)(a, k|s|a, k) = 0. In 
this way the derivative of energy with respect to mag- 
netization direction can be related to the exchange term 
in the Kohn-Sham equation rather than to the spin-orbit 
coupling term. Eqs. pO|) and (fTTj) are recovered after 
using [si, Sj\ = itijkSk to simplify ([A • s, s]). 

We now show that the Green's function expression we 
use to evaluate the ferromagnetic ISGE (Eq.([3])) corre- 
sponds to the current-induced change in Eq. (jlOp . We 
begin by mentioning that the application of an electric 
current can alter the magnetic anisotropy field, which 
leads to a current-induced torque on the magnetization. 
For an arbitrary orientation of the exchange field, the 
change is given by 
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<5H,, = ^ V S{{a, k| Ao(r)s|a, k»/k,a + ^T^ia, k| Ao(r)s|a, k)Sh,a 



r 



5*0 



(13) 



k,a 
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scattering rate r^^ in good conductors. However, as we 
show in the next section they may become quantitatively 
+... significant in disordered magnets hke the (III,Mn)V 
materials 1^ Admittedly, other corrections with the same 
parametric dependence on disorder strength could also 
be present - but the description of these would require a 
detailed characterization of the disorder potential and a 
more sophisticated transport theory. The effect we retain 
is analogous to the intrinsic contribution to the anoma- 
lous Hall effect,^ Substituting Eqs. dH]), ^ and Hr 
in Eq. we obtain 



FIG. 3: Spin response to a transverse magnetic field Bx in 
the presence of a current: perturbation theory to all orders in 
Bx- The quasiparticles (quasiholes) in these diagrams diag- 
onalize a Hamiltonian whose exchange field is pointing along 
an easy direction and B^ is by definition perpendicular to this 
easy direction. Provided that in Eq. (|10|l we take the exact 
eigenstates of the mean field Hamiltonian (within which the 
exchange field need not be pointing along an easy direction), 
all the diagrams of this figure are implicit in the diagram of 
Fig. IT}. In particular, the ferromagnetic ISGE captures the 
infiuence of currents on ferromagnetic resonance. 



where 



inter 
an 



5Han = SU" 



Ir V'[Ao(r)s]a,aVa,a ' E Tk,a 



k,a 

1 1 

iuj So 



[Ao(r)s]a,fcV6.a • E 



ek,& - ek,a +U! + 11] 



(17) 



Adopting the relaxation-time approximation, 6f reads 

(14) 



and for the change in the matrix elements we use 



S{{a, k| Aos|a, k)) = (a, k| AoS(5(|a, k)) + c.c (15) 



with 



<5(|a,k)) 



■E 



6,k) 



(6,k|v-E|a,k) 



Ek.a — ek,b 



+ {uj ^ -oj). 



(16) 

In Eq. (jl6p we have once again appealed to linear re- 
sponse theory and have used the fact that the electric 
field is uniform. 

Eqs. (fTil) and p6p highlight the two ways in which 
a current alters the magnetic anisotropy field. Eq. p4p 
captures the shift in the effective quasiparticle energies 
due to acceleration by an electric field, while Eq. pS]) 
describes the modification of the quasiparticle wavefunc- 
tions. As will become clear below the former is associ- 
ated with intraband contributions to the anisotropy field 
whereas the latter may be traced to the interband contri- 
butions. Interband contributions are often neglect ed^ii^ 
because they are parametrically smaller by a factor of 



In the expression for (JH^JJ*"' we have selected the co- 
efficient of exp(iwt) in the perturbation expansion, have 
neglected disorder scattering and have allowed for a small 
positive imaginary part in the frequency. 

Several remarks are pertinent in regards to our deriva- 
tion of the interband component. First, it should be 
noted that in the zero frequency limit the imaginary part 
of (5Hl^'Jj*°'^ gets cancelled by the diamagnetic contribution, 
in such a way that the anisotropy field induced by a DC 
current is finite and real. Second, it is instructive to 
elaborate on the real part of (JH^ 



rintcr . 
an 



SH 



inter 
an 



= 77— R-e [{ AqS) a, b^b, a] ' E(/k.a - fk.b)S{^b,a + w) 



k,a#& 



^ Im[(AoS)a,hVb a] ■E/k,a- 



2w 



(18) 



k.aj^b 



^b.a 



where LOh,a = Ck.fa — ek,a- /.From Eq. (|18p it is clear that 
(^HiJJj*""^ remains finite as w — > 0. When disorder is in- 
cluded in the above expressions, the contribution from 
the second line in Eq. scales as t^^ and thus is 

unimportant when the broadening of the energy bands 
due to impurity scattering is small compared to the en- 
ergy difference between states connected by interband 
transitions. In contrast, the third line scales as r'', and 
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therefore it supplies the bulk of the interband contribu- recall^ that 
tion in weakly disordered ferromagnets. 

Recognizing the fact that the integration of equal-band . . 

Green's functions gives rise to a factor of r, (JH!^"*'' ^ yields ^ ~ ^^'^ , ^ -^Y^ Ga(iw„, k)G6(iw„+it^, k). 

the intraband piece of Eq. ([2]) modulo a factor of Ao/5o. k ^^'^ ~ ~^ w„,k 

Similarly, (JHjjJ^"'' brings in the interband part of Eq. ^ (19) 
modulo a factor of Ao/5o; in order to verify this we In sum, we find 



^# = ^I^e E («,k|Ao(r)5>,k)(6,k|z;^|a,k) (G^,,Gi^,, - G«,G^,,) , (20) 

which agrees with the ISGE expression for the current-induced spin density (Eq. ^) except for an overall normal- 
ization factor (1/S'o) and the fact that the spin-operator is weighted by an spatially inhomogeneous magnitude of the 
exchange field. With the aim of making Eq. (j20p more manageable for first principles calculations, we will ignore the 
interband contribution as well as the G^G^ term; both omissions are justified in most metallic ferromagnetS;^ In 
this case Eq. (PD|) simplifies into 



an 



^E(«'k|^l«.k)(a,kK|a,k>%i^^rk,., (21) 



d5H._ 

k.a 



where we have re-inserted (a|Ao(r)s|a) = {a\d'Hso / dVt\a) . 
While approximate, Eq. (|2ip may provide a valid plat- 
form to explore current induced magnetization reversal 
in real gyrotropic ferromagnets with a single magnetic 
domain. In the next section we will describe in detail 
how a large i5ifan can produce a large reorientation of 
the magnetization. 

If the spatial dependence of Ao(r) is negligible (as it 
will be in the model studied in the next section) , Eq. (PD|) 
may be rewritten in a more compact way: 



Xs,E = 



Ao dEi 



(22) 



where xs,e is the spin-current susceptibility introduced 
in Eq. Eq. ([22]) proves that the ferromagnetic ISGE 
describes the change in the magnetic anisotropy field due 
to a current. In other words, ferromagnetic ISGE deter- 
mines how an electric current changes the location of the 
extrema in the micromagnetic energy functional. This is 
the central idea of this section. 

As a final sidenote, we point out that this section 
has concentrated on evaluating the change in magnetic 
anisotropy under a perturbation represented by v • A, 
where A is the electromagnetic vector potential. The 
anisotropy is evaluated by calculating the change in the 
expectation value of Aqs, thus leading to a rather stan- 
dard linear response function calculation. We could in 
the same way calculate the change in the transverse spin- 
spin response function due to an electric field as indicated 
in Fig.( [3]), in order to determine how small amplitude 
magnetic fluctuations are altered. If, however, we are 
interested only in uniform magnetization dynamics no 



additional information is obtained by doing this calcula- 
tion. The key point is that the response to a transverse 
field is already built in our expression for the equilib- 
rium anisotropy field (Eq. (HU])), to all orders in B^. In 
other words, the reference (unperturbed) macrostate to 
which we apply a current contains a magnetization that is 
"arbitrarily" misaligned with the exchange field. Hence, 
Eq. (|22p along with Eq. fTU]) offers a complete account of 
the nonequilibrium magnetic anisotropy of uniform mag- 
netic states in the presence of a transport current. 



III. CURRENT-DRIVEN MAGNETIZATION 
REVERSAL IN MONODOMAIN (GA,MN)AS 

Magnetoelectric phenomena in dilute magnetic 
semiconductors^ such as (Ga,Mn)As have attracted 
special attention because these materials are more com- 
patible with current microelectronics technology than 
metals. In addition, electric field control of magnetism 
has turned out to be more feasible in (Ga,Mn)As than 
in conventional dense-moment metallic ferromagnets 
because of their small magnetization, high carrier spin 
polarization, strong spin-orbit interactions, and carrier- 
mediated ferromagnetismi^i^i^ In particular, the recent 
experiment^ by Chernyshov et al. on (Ga,Mn)As wafers 
under compressive strain has demonstrated the ability of 
transport currents to reversibly assist the reorientation 
of magnetization in single-domain ferromagnets. As we 
demonstrate here this effect is dependent on having both 
spin-orbit interactions and broken inversion symmetry. 
In this section we compute the change in the magnetic 
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FIG. 4: Equilibrium anisotropy field (meV per spin) in 
(Ga,]V[n)As for = 0, and 9 £ (0,7r). The parameters used 
for this calculation were: Mn fraction x — 0.08, hole concen- 
tration p ~ O.lSnm^'', epT — 3, and axial strain eax = —0.5%. 
These anisotropy field results were evaluated using the model 
explained in the text. 
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FIG. 5: Equilibrium anisotropy field (meV per spin) in 
(Ga,Mn)As for 9 — tt/2 and £ (0, vr). The parameters 
are: Mn fraction x — 0.08, hole concentration p ~ 0.15nm~^, 
epT = 3, and axial strain tax = —0.5%. These results were 
evaluated using the model explained in the text. Due to 
strain, the in-plane anisotropy is notably weaker than the 
out-of-plane anisotropy represented in the previous figure. 



FIG. 6: Change in the magnetic anisotropy field of 
(Ga,Mn)As (in meV per spin) due to the inverse spin-galvanic 
effect, for an electric field of ImV/nm along [010]. The 
parameters are: Mn fraction x — 0.08, hole concentration 
~ 0.25nm~^, epr — 2, and axial strain tax ~ —1% We com- 
pare between interband and intraband contributions: in con- 
trast to the case of good metals, the interband contributions 
are not negligible in (Ga,Mn)As. For the present case, had 
we neglected the interband contribution the minimum electric 
field needed to reorient the magnetization by 90° would be off 
by approximately 20 %. The sum of interband and intraband 
pieces gives rise to a smooth curve that reflects the Dressel- 
haus symmetry of the axial strain. Reversing the sign of the 
axial strain (i.e. malting it tensile) leads to a sign reversal of 



ture and strain has been suecessfuUy explained^i^i^ by 
combining (i) a mean-field theory of the exchange cou- 
phng between locahzed Mn moments and valence band 
carriers with (ii) a phenomenological four or six band en- 
velope function model in which the valence band holes 
are characterized by Luttinger, spin-orbit splitting and 
strain-energy parameters. The results presented below 
predict the rate at which these fields change with exter- 
nal electric field. 

In line with this we adopt the following Hamiltonian 
for Gai_3;Mn2:As: 



H = Hkl 



n. 



St ram 



S- A. 



(23) 



anisotropy due to an electric current for a realistic model 
of (Ga,Mn)As. Our calculation is directly relevant 
to the experiment of Chernyshov et ai. Our results 
corroborate their interpretation of the data and predict 
the possibility of all-electric magnetization switching in 
(Ga,Mn)As. Our analysis is limited to zero temperature 
and neglects the shape anisotropy, which for typical Mn 
doping concentrations is 10-100 times weaker than in 
conventional fcrromagncts. 

The dependence of the magnetic anisotropy of 
(Ga,Mn)As on doping, external electric fields, tempera- 



Hkl is the 4-band Kohn-Luttinger Hamiltonians^ with 
Luttinger parameters 71 = 6.98, 72 = 2.1 and 73 = 2.9. 
S is the spin operator projected onto the 3=3/2 total 
angular momentum subspace at the top of the valence 
band. A = AqCi = Jp^S Nun^ is the exchange field, 
f2 denotes the orientation of the local moments, Jpd = 
55 mcV nm is the p-d exchange coupling parameter, S = 
5/2 is the spin of the Mn ions, and A^Mn = 4a;/a^ is the 
Mn concentration (a = 0.565nm is the lattice constant of 
GaAs). This four-band model is expected to be adequate 
for small and intermediate Mn doping strengths. Hstrain 
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FIG. 7: Reorientation of the magnetization due to an electric 
current. An initial magnetization along [100] can be rotated 
(assisted by damping) into [010] by applying a sufficiently 
strong electric field with a nonzero projection along the [010] 
direction (a current along [100] would not destabilize the [100] 
easy axis). For the parameters of this figure {x = 0.08, p ~ 
0.15nm~'^, epT = 3, tax = —0.5%) the critical electric field is 
~ 5mV/nm, which corresponds roughly to a critical current 
density of 5 x lO^A/cm^. 



FIG. 8: Dependence of the critical electric field (at which the 
magnetization gets reoriented by 90°) on (compressive) axial 
strain. The critical current is (roughly) inversely proportional 
to tax. The reason behind this relationship is that the equi- 
librium, azimuthal anisotropy is largely indiferent to tax. For 
X = 0.04 and tax = —2% we find Ec ~ 0.25mV/nm, which 
corresponds to a critical current on the order of lO^A/cm^. 
These results are for a (Ga,Mn)As model with carrier density 
p ~ 0.15nm~'' and eft = 3. 



is the strain Haniiltonian^ i"^^'^^ given by 



strain 



= -b 



c.p. 



+ Ci [Jx {eyy - £zz) + c.p] 



(24) 



where J is the total angular momentum (J = 3S by 
the Wigner-Eckart theorem), ti^i are diagonal elements 
of the stress tensor, b = —1.7 eV is the axial deforma- 
tion potential and the parameter C4 = 5 eV A captures 
the strain-induced linear in k spin-splitting of the valence 
bands in paramagnetic GaAs. In Eq. the notation 
c.p. stands for cyclic permutations and ex,x = ^y,y 7^ ^z,z 
for [001] growth lattice-matching strains. The term pro- 
portional to C4 is crucial for the occurrence of the fer- 
romagnetic ISGE because it breaks inversion symmetry 
(we are neglecting the intrinsic lack of inversion symme- 
try of the zinc-blende structure, which is relatively in- 
consequential), and it introduces chirality. (A bulk, un- 
strained zinc-blende crystal is not gyrotropic because it 
corresponds to the Td symmetry point group.) Eq. (j24[) 
may be simplified to 



"^strain ^^a 



Jz ) + C4£ax {Jyky — Jxkx) ■ 



(25) 

where eax = Czz — ^xx is the purely axial strain compo- 
nent. In this paper we take tax < (compressive strain), 
which applies when (Ga,Mn)As is grown on top of a GaAs 
substrate. 

Using Eqs. HU]) , (EH) and ^ we evaluate the 
magnetic anisotropy field both with and without elec- 
tric current; the results are highlighted in Figs. ([3])- ([5]). 



Figs. Q and ([5]) correspond to electrical equilibrium and 
illustrate Hg = -1/S'o I]k,a(^'^k,Q/96')/k,a for = and 
= -l/S'oX]k,a(^'^k,a/90)/k,a for 9 = tt/2, respec- 
tively. The extrema of the micromagnetic energy func- 
tional are characterized by = Hg = and by inspec- 
tion we locate them at 9 = and {9,(j)) = (7r/2,n7r/4) 
where n ~ 0, 1,2.... For our parameters (see figure cap- 
tions) the energy minima that define metastable mag- 
netic configurations are found at {9,(j)) = (7r/2, n7r/2). 
That is to say, the easy directions correspond to [100], 
[010], [100] and [010], which are contained in the plane of 
the (Ga,Mn)As wafer. For later reference, we consider an 
initial condition in which the magnetization is pointing 
along [100]. If a small static perturbation tilts it towards 
[110], the negative anisotropy field {H^ < for > 0) 
creates a torque that will, in conjunction with damping)^ 
turn the magnetization back to [100]. 

Fig. ([6]) illustrates how an electric current along [010] 
alters the azimuthal anisotropy field^ for fixed 9 = tt/2. 
The cosine-like shape is consistent with the Dresselhaus 
symmetry of the C4 term in the strain Hamiltonian. 
If the system had a perfect Dresselhaus symmetry the 
change in the micromagnetic energy functional under an 
electric current j would read 



SEgs « Ciei,^{^yjy - il.xjx), 



(26) 



which means that a current along [010] ([100]) would 
tilt the steady-state magnetization direction along [010] 
([100]). Using = sin0cos(/) and fiy = sin sine/) it fol- 
lows that SH^ oc jy cos 4'+jx sin 0, and hence a cosine- like 
dependence in (f> is indeed expected for a current along 
[010] . We have verified that a current along x gives rise to 
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a sine-like dependence with the appropriate sign. Nev- 
ertheless, Eq. (|26p is not exact because the magneti- 
zation vector introduces another preferred direction; for 
instance, we find that an electric field pointing along z 
(i.e. [OOf]) can also alter the steady-state spin orienta- 
tion. This effect, which vanishes in the paramagnetic 
limit, highlights one instance in which the ferromagnetic 
and paramagnetic ISGEs differ. Another attribute of Fig. 
([61) is that it determines the quantitative importance of 
interband contributions to the current-induced spin den- 
sity in (Ga,Mn)As. Although normally neglected, in- 
terband transitions become quantitatively significant in 
strongly disordered ferromagnets. fn particular, inter- 
band and intraband contributions are largely indistin- 
guishable in ferromagnets with Aqt < 1. We note par- 
enthetically that neither intraband nor interband contri- 
butions display the smooth sinusoidal shape portrayed 
by their sum. In addition, we remark that reversing the 
sign of the axial strain (i.e. making it tensile) leads to 
a sign reversal of bH^ without substantial changes in its 
magnitude?^ 

Fig. ([7]) demonstrates that a sufficiently strong current 
is able to rotate the magnetization by 90° or f 80°. We 
explain this property by considering the case in which 
the equilibrium magnetization is pointing along [100]. If 
a small current is applied along [010], then [100] is no 
longer an extremum of the micromagnetic energy func- 
tional (because = 0) oc i?y ^ 0). The modified easy 
direction remains in the neighborhood of [100] since the 
restoring torque {H^ < 0) again crosses zero at (j> > 0. 
Once the applied electric field exceeds a critical value 
{Ec — 5.5mV/nm in the present figure) the < re- 
gion near [100] disappears completely and hence assisted 
by damping the magnetization eventually points along 
[010]. In other words, at (and above) the critical switch- 
ing field the energy minimum that is nearest to [100] is 
located at [010] (note that this direction remains stable 
when the current flows along [010]). Once the magne- 
tization is aligned with [010], an equally strong electric 
current in the [100] direction will rotate it towards [100]. 
In this fashion it is possible to switch the direction of 
magnetization by 180° solely by application of transport 
currents. 

The procedure sketched above accomplishes magneti- 
zation switching by application of two perpendicular cur- 
rent pulses, each of which forces a 90° rotation. Yet, it is 
also possible to achieve the [100] [100] switching with 
a single unidirectional pulse, provided the electric field 
along [100] is ramped up sufficiently (i?c,2 — 20mV/nm 
for the parameters of the present figure). In order to 
understand this, recall that j||x ^ <5Han|| — x. Conse- 
quently, for a strong electric current [100] is the only 
easy direction ([100] becomes a hard direction). The 
inequivalence between [100] and [100] does not violate 
any symmetry principles'^ in effect, an electric current 
breaks time reversal symmetry and can thus connect 
time-reversed magnetic states. 

Using p ~ 10^'^Ocm as the typical resistivity for 



(Ga,Mn)As samples we deduce that E = ImV/nm corre- 
sponds approximately to a current density of lO^A/cm^, 
hence the critical switching current is on the order of 
10® — lO^A/cm^. It is plausible that a detailed explo- 
ration of the parameter space comprised by the Mn con- 
centration X, the hole density p and the axial strain tax 
will enable lower critical currents, thereby diminishing 
the importance of the Joule heating. As a word of cau- 
tion, we note that the 4-band model employed here typ- 
ically overestimates the effect of spin-orbit interactions, 
thus potentially leading to an underestimate of these crit- 
ical currents. There is in addition some uncertainty asso- 
ciated with the use of a life-time approximation for Bloch 
state quasiparticles in these strongly disordered metallic 
conducting ferromagnets. 

Overall, the magnitude of the critical switching cur- 
rent depends on (a) the size of the equilibrium anisotropy 
barrier, (b) the extent to which inversion symmetry is 
broken and (c) the strength of spin-orbit interaction. 
In (Ga,Mn)As the first two factors are tunable. On 
one hand, (a) may be optimized by choosing appropri- 
ate doping concentrations: in general lower Mn den- 
sity is beneficial (Fig. as it reduces the equilib- 
rium anisotropy without significantly affecting the mag- 
nitude of ISGE. However, for very low Mn concentrations 
a metal-insulator transition is impending, which ham- 
pers ISGE. On the other hand, (b) may be modified via 
strain engineering: as shown in Fig. ([8]), the critical cur- 
rent is (roughly) inversely proportional to the strength of 
the uniaxial strain that breaks inversion symmetry. The 
inverse proportionality may be understood on the basis 
of Eq. (|26|) combined with the fact that the equilibrium 
anisotropy does not change to first order in eax (because 
k-linear terms vanish after summing over all momenta). 

IV. SUMMARY AND CONCLUSIONS 

In this work we have presented a theory of the current- 
induced spin torques in uniform ferromagnets. The 
torques can be viewed as due to a difference between 
the magnetic anisotropy energy of a ferromagnet which 
carries no current and the magnetic anisotropy of a ferro- 
magnet in the transport steady state, which give rise to 
a corresponding change in anisotropy effective magnetic 
fields. When the transport steady state is described us- 
ing a relaxation time approximation, the current-induced 
contribution to the anisotropy field of a strongly metallic 
ferromagnet is given in energy units by 

<5Han = ^ V[Ao(r)s],,,v„,„ • E |^ r^.a- (27) 

k,a 

where [Ao(r)s]a,a is the spin-density weighted average of 
the exchange splitting of a particular state. We refer to 
the existence of this current-induced anisotropy field as 
the ferromagnetic inverse spin-galvanic effect. 

In bulk materials this current induced field is non-zero 
only in gyrotropic ferromagnets, i.e. only in ferromag- 
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nets that have broken inversion symmetry and are chi- 
raL Although uniform ferromagnetism may appear to be 
incompatible with broken inversion symmetry because of 
the the Dzyaloshinskii-Moriya interaction, the equilib- 
rium magnetic anisotropy is often strong enough (or at 
least can be engineered so that it is strong enough) to 
prevent the formation of spiral magnetic states. 

As an illustration of our theory, we have estimated 
current induced torques in uniform (Ga,Mn)As, which 
is not gyrotropic when it has pseudo-cubic symmetry 
but becomes gyrotropic when strained. Since substrate- 
dependent strains are present in all (Ga,Mn)As thin 
films, the strength of the ferromagnetic ISGE is expected 
to be strongly sample-dependent. We have concluded 
that it should a priori be feasible to design (Ga,Mn)As 
samples in which it is possible to switch the magneti- 
zation purely by electrical means. For typical sample 
parameters the necessary switching currents are on the 
order of 10^ — lO^A/cm^, but the value may be tuned by 
adjusting the doping concentration and the axial strain. 
At these critical currents the Joule heating is not negli- 
gible; however, it is possible that further studies explor- 
ing the entire parameter space of Mn concentration, hole 
density, and the axial strain will identify circumstances 
under which the critical currents are smaller. 



Another possible avenue for further research consists 
of evaluating the anisotropy fields which can be gen- 
erated by electrical currents in strain engineered sam- 
ples of appropriate technologically useful ferromagnets. 
Since we are not aware of room-temperature transition 
metal ferromagnets that are gyrotropic,— we propose ar- 
ranging a room-temperature, non-gyrotropic ferromag- 
net (e.g. permalloy) in contact with a non-magnetic, 
gyrotropic material (e.g. strained GaAs). In these arti- 
ficial heterostructures room-temperature magnetism and 
gyrotropic symmetry would coexist by virtue of the prox- 
imity effect. 

Finally, effects similar to those studied in this work 
would allow transport currents to change spiral states, 
and possibly to induce or remove them. 
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